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Abstract: Diffuse reflectance spectroscopy (DRS) is a well-established 
method to quantitatively distinguish between benign and cancerous tissue 
for tumor margin assessment. Current multipixel DRS margin assessment 
tools are bulky fiber-based probes that have limited scalability. Reported 
herein is a new approach to multipixel DRS probe design, which utilizes 
direct detection of the DRS signal by using optimized custom 
photodetectors in direct contact with the tissue. This first fiberless DRS 
imaging system for tumor margin assessment consists of a 4 x 4 array of 
annular silicon photodetectors and a constrained free-space light delivery 
tube optimized to deliver light across a 256 mm 2 imaging area. This system 
has 4.5 mm spatial resolution. The signal-to-noise ratio measured for normal 
and malignant breast tissue-mimicking phantoms was 35 dB to 45 dB for X 
= 470 nm to 600 nm. 
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1. Introduction 

Diffuse reflectance from turbid media such as biological tissue is sensitive to the absorption 
and scattering of the constituents of these media. In the case of biological tissue, diffuse 
reflectance spectroscopy (DRS) can be used to look for spectral signatures of molecules such 
as oxy -hemoglobin, deoxy-hemoglobin, collagen, and beta-carotene [1,2], for example. DRS 
performed with well-defined illumination and collection geometries and coupled with 
physically based light transport models can be used to quantitatively extract tissue optical 
properties (absorption and scattering) [3-5]. This information can then be used to develop 
appropriate decision tree models to discern between healthy and cancerous tissue [2,6]. 

One application for which DRS has been investigated is intra-operative tumor margin 
assessment for breast conserving surgery (BCS) [2,6,7]. Over 200,000 women are diagnosed 
with breast cancer annually in the United States, and 50%-75% of these women are eligible 
for BCS [8]. Of the women undergoing BCS, 20%-70% return for a re-excision surgery due 
to incompletely removed tumor in the first surgery [9-11]. The goal of DRS imaging in this 
context is to reduce the re-excision rate. 

Our group has developed a bench-top tissue imaging system, based on quantitative DRS, 
for intra-operative tumor margin assessment. Briefly, the bench-top system consists of a 
xenon arc lamp, an imaging spectrograph, and a cooled CCD camera to image the DRS 
spectra collected by an 8-channel fiber optic probe [12]. The 8-channel fiber optic probe 
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consists of 8 individual fiber probes, each composed of illumination fibers surrounded by 
collection fibers. An inverse Monte-Carlo model, developed by Palmer et al. [4], is used to 
extract tissue optical properties (absorption coefficient and reduced scattering coefficient) 
using the DRS spectra collected by the bench-top system. The system surveys the target 
region of interest over a wavelength range of 450-600 nm and a sensing depth of 2 mm (the 
criterion for clear margins) to provide quantitative maps of oxy-hemoglobin, deoxy- 
hemoglobin, beta-carotene and scattering. A classification model was developed using a 
combination of these endpoints to detect positive tumor margins with a sensitivity of 79.4% 
(i.e., percentage of correctly detected positive margins), and specificity of 66.7% (i.e., 
percentage of correctly detected negative margins) from a 48 patient cohort, demonstrating its 
potential to significantly reduce re-excision surgeries [2]. However, this system costs $50,000, 
has a footprint of 2 m x 1 m, and is only capable of coarse 8-pixel imaging [12,13]. In 
addition, the ability to rapidly assess large margin areas is key in the context of breast cancer 
imaging, since excised breast tissue areas can be as large as 20 cm 2 . Therefore, there is a need 
for a low-cost, fast, compact, portable spectral imaging system with multi-pixel, large-area 
probes for breast tumor margin assessment inside the surgical suite, where space and time are 
at a premium. 

In an effort to reduce the complexity and footprint of this bench-top system, a more 
compact system was developed where DRS spectra were sampled using a xenon lamp 
followed by eight discrete wavelength filters. It was demonstrated that optical properties could 
be extracted with low errors using the inverse Monte Carlo model, with eight wavelengths of 
light as opposed to 81 wavelengths used in the bench-top system [14]. In addition, an array of 
mechanically drilled commercial silicon photodiodes was developed to achieve Si 
photodiode-based collection, and fibers were epoxied to the drilled photodiode apertures to 
achieve fiber-based illumination in the DRS probe [13]. This system obviated the need to use 
an imaging spectrograph and a 2D CCD for obtaining spectral information; components that 
significantly added to the complexity and foot-print of the bench-top system. The drilled 
photodiodes with epoxied fibers replaced previous illumination/collection fiber probes, 
however pixel spacing was still limited by commercial photodiode packaging and this system 
had limited scalability and resolution. 

Here, we report the design, optimization, fabrication, and characterization of a multipixel 
DRS imaging system using a custom photodetector array and constrained free space 
illumination. This system enables fiberless illumination of tissue sites and collection of DRS 
signals by the use of an array of custom annular silicon photodetectors (PDs), which directly 
contact the tissue. Illumination from a filtered xenon lamp is incident upon an aperture in the 
center of each of the annular PDs. After the signal enters tissue and is diffusely reflected, the 
DRS signal is collected by the annular PD surrounding the aperture. This illumination- 
collection geometry has been shown to have the highest SNR for contact -based DRS [14,15]. 
While annular silicon PDs with back illumination have been used in applications such as pulse 
oximetry [16], these systems utilize single PDs and have two -wavelength interrogation. This 
report utilizes an array of sixteen annular PD pixels, arranged in a 4 x 4 format and eight 
discrete wavelengths to obtain DRS spectral information from thick tissues. 

2. System description 

The schematic of the overall system is shown in Fig. 1(a). Figure 1(b) shows a photograph of 
the custom compact spectral imaging system, and Fig. 1(c) shows the contact imaging probe. 
The system reported herein consisted of a xenon arc lamp (MAX-302, Asahi Spectra) with 
eight discrete wavelength filters as the illumination source, a sixteen pixel imaging array 
composed of customized annular silicon photodiodes, a quartz light delivery tube that enables 
free space illumination of tissue sites, a custom eight-channel integrating transimpedance 
amplifier (ITIA) circuit for photodiode readout, and a custom switch box connected to the 
eight channel ITIA circuit that enables per-pixel amplification and readout of all sixteen pixels 
in the array. Each of these system components is controlled by a Lab VIEW program running 
on a laptop computer. 
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Spectral information was collected by making successive diffuse reflectance 
measurements with each of the eight wavelength band pass filters with center wavelengths at 
470 nm, 480 nm, 490 nm, 500 nm, 510 nm, 560 nm, 580 nm, and 600 nm [17]. Each filter was 
mounted on an eight-slot filter wheel inside the xenon arc lamp and had a AX — 10 nm full 
width at half max. Light from the xenon lamp passed through one filter at a time, and exited 
the lamp through a 5.2 mm diameter, 0.57 numerical aperture (N.A.) fiber bundle. The fiber 
bundle was connected to one end of the quartz light delivery tube and the custom imager was 
mounted on the distal end of the quartz tube, as shown in Fig. 1(a) and Fig. 1(c). The imager 
consisted of sixteen annular silicon photodiodes with 2.5 mm outer diameter (O.D.), and 0.75 
mm inner diameter (I.D.) central apertures. The center-to-center spacing of the array pixels 
was 4.5 mm. The O.D. and I.D. of the PDs were designed to achieve a sensing depth of 0.4-2 
mm in breast tissues [18]. Light traveled from the fiber bundle, through free space in the 
quartz tube, through the PD apertures. The specimen under investigation was placed in contact 
with the imager. Light traveling through the PD apertures illuminated the tissue sites, and 
diffusely reflected light from the specimen was incident upon the detection area of the annular 
PDs. This diffusely reflected light generated carriers in the photodiodes, read as photocurrent 
by the readout circuit. 




(b) (c) 



Fig. 1 . (a) Schematic of the DRS imaging system and cross-section view of a single PD in the 
array (b) Photograph of DRS imaging system (c) Photograph of the customized DRS imaging 
probe with a fiberless quartz light delivery tube and an annular custom silicon photodiode 
array. 

3. System development and optimization 

3.1. Custom annular photodetector array fabrication and characterization 

Fabrication of the custom-imaging array began with an n-type epitaxial Si wafer: 7 urn «-Si 
(epitaxial, P doped, 5-9 Q-cm) on 600 urn «-Si (substrate, Sb doped, 0.025-0.045 Q-cm). The 
Sb-doped substrate side of the wafer was coated with a phosphorous doped spin-on-glass 
(Emulsitone Phosphorosilica film 5 x 10 20 cm -3 ) and annealed to form an n + back contact 
layer. A 150 nm thick Si0 2 layer was grown using thermal oxidation on the epilayer and was 
patterned to define the active area of each PD. A boron doped spin-on-glass (Emulsitone 
Borosilicafilm 5 x 10 19 cm -3 ) was spin coated on the front-side epitaxial layer of the wafer 
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and annealed to form the pn junctions. Since the diffusion coefficient of boron is lower in 
Si0 2 compared to Si [19], the patterned oxide layer served as a diffusion mask to dope only 
the active areas of the PDs. The top dual-ring aluminum ohmic contact was patterned using 
UV lithography and an Al/Pt stack was deposited as the blanket ohmic back contact. A 0.75 
mm diameter 600 um deep aperture was patterned in each PD using UV lithography and 
etched using deep reactive-ion etching (Bosch process, SPTS Pegasus). Finally, a 52 nm thick 
PECVD silicon nitride (SiN) layer was deposited on the PDs as an anti-reflection (AR) 
coating. The SiN layer also served to passivate the front surface of the PDs, thereby reducing 
the front surface recombination velocity [20]. The fabricated imaging array was wire-bonded 
to a patterned gold-coated FR4 board to enable PD current readout. 

Characterization of the imaging array was performed using the fiber-coupled Asahi lamp 
with 8 wavelength filters as an illumination source to measure surface normal spectral 
response of the fabricated PDs. Light from a fiber bundle was incident on each PD in the 
surface normal configuration, and the resulting photocurrent in each PD was measured using a 
source-measurement unit (Keithley SMU-4200). Power at the output of the fiber bundle was 
measured using a calibrated Si photodetector head (Thorlabs S120A) connected to an optical 
power meter (Thorlabs PM100) to calculate the PD spectral response. A 3D simulation was 
performed in Silvaco Atlas® to theoretically estimate the surface normal spectral response 
and junction depth of the fabricated PDs. The Silvaco model used the dopant concentrations 
and epilayer thickness of the starting wafer. A complementary error-function profile was 
assumed for boron diffusion with boron surface concentration equal to 5 x 10 19 cm -3 . Figure 
2(a) shows the average measured surface normal spectral response for 15 PDs in the array 
shown in Fig. 2(b). One pixel in the array had an unusually low responsivity of 1 mA/W to 3 
mA/W, and was not included in the average. The simulated spectral response for a junction 
depth of 700 nm agrees well with the measured responsivity. The average responsivity for the 
15 PDs was measured to be 0.13 AAV to 0.3 AAV for X = 400 nm to 600 nm. The dark current 
density of the fabricated PDs was measured to be 14.56 pA/mm 2 to 44.8 pA/mm 2 . 

■ Measured spectral response 



» Simulated spectral response: junction depth = 700 nm 




1 1 1 1 1 1 1 1 1 1 1 1 

350 400 450 500 550 600 650 

Wavelength (nm) 

(a) (b) 

Fig. 2. (a) Measured and predicted spectral response of fabricated 4x4 array of silicon PDs. 
Error bars are ±la. (b) Photograph of fabricated and packaged array (inset) Photomicrograph of 
one pixel in the annular PD array. 

3.2. Light delivery design and characterization 

Light was delivered to individual test sites through central apertures in the fabricated PDs to 
maximize the SNR of the resulting diffuse reflectance signal [15]. While our previously 



#176730 - $15.00 USD Received 27 Sep 2012; revised 7 Nov 2012; accepted 7 Nov 2012; published 8 Nov 2012 
(Q2012OSA 1 December 2012 /Vol. 3, No. 12 / BIOMEDICAL OPTICS EXPRESS 3215 



reported diffuse reflectance imaging systems have utilized individual fibers or fiber bundles to 
illuminate each tissue site, these systems have been limited to less than 10 illumination sites 
per image, and pixel spacings greater than 7 mm [12,13]. As pixel to pixel spacing decreases, 
and in addition, imaging area increases, illuminating each tissue site with a fiber becomes 
increasingly unwieldy. One way towards a more scalable light delivery system without fiber- 
based illumination is by using a light delivery tube coupled to a single source fiber bundle, to 
deliver light to all of the pixels of the array. 

The system reported herein uses a metalized quartz tube, optimized in length and diameter, 
to ensure sufficient power throughput through PD apertures and uniformity across all pixels. 
This optimized tube replaced a suboptimal Al tube used in previous experiments [21,22]. The 
diameter of the fabricated array is 22.6 mm, so a 27 mm I.D/30 mm O.D. quartz tube was 
selected, in order to accommodate the array in the center of the tube for illumination. The 
length and diameter of the tube were simulated to utilize rays reflected from the metalized 
tube to enhance light throughput and uniformity across the array. Simulations were performed 
using a commercially available ray tracing software, ZEMAX®, to estimate light throughput 
and uniformity. The ZEMAX® model simulated the NA. (0.57) of the source fiber bundle 
and the intensity distribution across the source fiber bundle [23]. A tube diameter larger than 
30 mm O.D. increases the length of the tube required to utilize reflections from the metalized 
walls of the tube without any improvement in light throughput. Therefore, in order to limit the 
size of the probe, a tube diameter of 27 mm I.D./30 mm O.D. was chosen. Optimization of the 
tube length, for the chosen tube diameter, is shown in Fig. 3(a) for X = 470 nm and 600 nm. It 
should be noted that corner pixel illumination is a challenge due to the thickness of the PD 
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Fig. 3. (a) Power throughput through PD apertures using a reflective tube for light delivery, 
simulated in ZEMAX® (b) Average measured power throughput for 4 center pixels, 4 comer 
pixels and 8 side pixels using a 70 mm long reflective tube (c) Average measured power 
throughput for 4 center pixels, 4 corner pixels and 8 side pixels using a 70 mm long absorbing 
tube. Error bars are +lo, and the spectral content reflects that of the source. 
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array and their distance from the center of the fiber bundle. A tube length of 70 mm is 
optimum for corner pixel throughput and uniformity between corner and center pixel light 
throughput. The simulated corner-to-center pixel light throughput ratio was 0.78 (i.e., 78% 
uniformity) for a 70 mm tube length. 

The quartz tube was coated with Ag (2500 A)/Au (2500 A) using e-beam evaporation, due 
to high reflectivity of Ag for X = 470 nm to 600 nm [24]. The measured optical power through 
the PD apertures is shown in Fig. 3(b). The measured corner-to-center pixel uniformity was 
73.5%, and was within 5% of the value predicted by ZEMAX®. 

Power throughput through the PD apertures was also measured for a quartz tube identical 
in dimensions to the reflective tube, but coated on the inside with absorbing black paint. 
While using an absorbing tube reduces light throughput by 30%, as shown in Fig. 3(c), the 
angular distribution of rays emerging from the PD apertures is smaller for the absorbing tube 
compared to the reflective tube (due to the absence of reflected rays in the absorbing tube). 
The angular distribution of rays has an impact on pixel-to-pixel optical crosstalk, with 
increasing angular distribution leading to higher levels of crosstalk, which has also been 
studied herein. 

3.3. Custom interface circuit design and development 

Eight parallel Integrating Transimpedance Amplifiers (ITIAs) were built using commercial 
parts [25]. A digitally controlled mechanical relay switch box was used to connect the 8 
amplifiers to the 16 pixels, 8 at a time. A low cost microcontroller [26] was used to control the 
amplifier operation and interface to a measurement computer and data collection system. The 
interface circuit was programmable for different integration times, and this system achieved 
an SNR of better than 55 dB for currents greater than 4 nA with 100 ms integration time. The 
entire resulting interface system, including ITIAs, switch box, and data acquisition, occupies a 
volume that is significantly smaller than the volume occupied by the xenon lamp. 

4. Experimental system characterization 

Experimental characterization of the system was conducted to quantify system signal-to-noise 
ratio (SNR) and optical crosstalk in tissue. Optical crosstalk in tissue is the diffuse reflectance 
signal generated in a pixel due to tissue illumination from neighboring pixel apertures. System 
SNR and crosstalk were measured for two breast tissue mimicking phantoms. The phantoms 
were a mixture of methemoglobin (absorber) and polystyrene spheres (scatterers) dissolved in 
deionized water. The phantoms were prepared with known optical properties similar to normal 
and malignant breast tissue. The normal tissue phantom is characterized by low absorption 
and high scattering properties (averaged (x a = 0.56 cm -1 and \i s ' = 11.9 cm -1 ), whereas the 
simulated malignant tissue has relatively higher absorption and lower scattering (averaged 
optical properties (x a = 5.42 cm -1 and |x s ' = 3.43 cm -1 ). Because the goal is to study the effect 
of optical crosstalk in tissue, the worst-case and best-case scenarios for optical crosstalk in 
clinically relevant tissue were simulated based on previously measured breast tissue data [27]. 
The assumption is that any mean or median normal and malignant tissue types would not 

Table 1. Optical properties of tissue phantoms 



Benign breast tissue Cancerous breast tissue 
Wavelength phantom [27] phantom [27] 



(nm) 


ft,, (cm ') 


fi s '(cm ) 


fi a (cm ) 


fi s '(cm ) 


470 


0.6646 


11.1447 


6.5732 


3.1904 


480 


0.6465 


10.9422 


6.3943 


3.1324 


490 


0.6603 


10.8254 


6.5300 


3.0990 


500 


0.6644 


10.7530 


6.5708 


3.0783 


510 


0.6307 


10.6877 


6.2380 


3.0596 


560 


0.4254 


10.0198 


4.2068 


2.8684 


580 


0.4207 


9.7991 


4.1611 


2.8052 


600 
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exhibit as much optical crosstalk, and it is important to characterize and model the crosstalk in 
normal and malignant tissues to be able to ultimately classify them. The wavelength- 
dependent optical properties of the tissue -simulating phantoms are listed in Table 1. 

4.1. SNR 

The system SNR was calculated using 30 repeated measures at each wavelength of diffuse 
reflectance signal photocurrent generated in each PD when the imaging array was placed in 
contact with the tissue phantom. The SNR in dB is defined as 20log [0 (I mean /(j), where I mean is 
the mean of 30 repeated measures, and a is the corresponding standard deviation. 
Unfortunately, in the course of imaging array characterization, two pixels in the imaging array 
became unresponsive due to disconnected wire-bonds, but these poorly responding pixels are 
simple packaging and yield problems that will be addressed in future work. Figure 4(a) and 
Fig. 4(b) show the average SNR for 13 out of 16 pixels in the imaging array (excluding the 
two unresponsive pixels, and one pixel with low responsivity), measured from the normal 
tissue phantom tested with the reflective tube for light delivery, and for a malignant tissue 
phantom tested with the absorbing tube for light delivery. 
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Fig. 4. (a) Calculated SNR for normal tissue phantom tested with a reflective tube for light 
delivery using 10 ms integration time in ITIA read-out circuit; (b) Calculated SNR for 
malignant tissue phantom tested with an absorbing tube for light delivery using 100 ms 
integration time in ITIA readout circuit. Error bars are +la for the mean SNR across 13 pixels. 

It is observed from Fig. 4(a) and Fig. 4(b) that higher SNR is achieved for the malignant 
tissue phantom measured using the absorbing tube in spite of lower signal levels compared to 
normal tissue phantom (less absorbing) measured using the reflective tube (higher 
throughput). These higher signal levels from the reflective tube dictated shorter integration 
times, to avoid saturating the TIA. The higher SNR from the absorbing tube is likely due to 
lower o as a result of averaging of the lamp noise with the longer integration time (100 ms). 

4.2. System crosstalk measurement and analysis 

In the design of fast multipixel DRS probes, it is important to understand the trade-off 
between resolution (pixel to pixel spacing) and crosstalk. Optical crosstalk in tissue refers to 
the photocurrent generated in a PD pixel due to illumination of tissue from neighboring PD 
apertures. A smaller pixel to pixel spacing (i.e., higher resolution) results in higher crosstalk. 
As the pixel to pixel spacing (i.e., PD spacing) in the imaging array decreases, more photons 
can travel from through tissue from neighboring PD apertures. This may reduce contrast 
between two neighboring PD pixels because the photons detected by a PD can no longer be 
correlated to the illumination aperture of the same PD. Hence, optical crosstalk in tissue is an 
important system parameter for multipixel DRS imaging systems. In this report, crosstalk is 
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defined as the percentage of total diffuse reflectance signal generated at a pixel due to photons 
traveling through tissue from neighboring illuminated PD apertures. In previously reported 
DRS systems, crosstalk comprised less than 10% of the total measured diffuse reflectance 
signal at each pixel. Crosstalk can be reduced by increasing center-to-center pixel spacing, as 
in these previously reported systems, which had a center-to-center pixel spacing greater than 
or equal to 7 mm [12,13], in contrast to the 4.5 mm center-to-center pixel spacing reported 
herein. 

Optical crosstalk in tissue was simulated for the probe geometry reported herein, using a 
scalable forward Monte-Carlo model previously described [28]. Photon launch trajectories 
were first simulated by modeling the illumination geometry in ZEMAX®. This provided the 
ray angles emerging through each aperture in the PD array. The optical power through each 
aperture in the array was measured and used to model the relative intensity of each aperture in 
the array. The ray angles, along with the relative aperture intensities, were used in a set of 16 
different baseline simulations, corresponding to each aperture in the array. The boundary 
conditions were assumed to be semi-infinite, with a 52 nm thick layer SiN coating over 
silicon, whose reflectivity was calculated using the method described in Born and Wolf [29]. 
All simulations were performed using optical properties of SiN calculated for X = 500 nm 
[30]. The signal generated at a PD due to photons traveling through its own illumination 
aperture was determined by calculating the collected fraction hitting each detector when its 
own aperture was illuminated, while the signal from the illumination of neighboring apertures 
was calculated by summing the contribution from the remaining 15 apertures. 

Figure 5(a) and Fig. 5(b) show the simulated optical crosstalk for a center pixel of the 
imaging array using a reflective tube, and an absorbing tube, respectively. The contour lines 
represent log ^crosstalk), and red and green data points represent the theoretical optical 
properties of malignant and normal tissue mimicking phantoms (listed in Table 1) at the eight 
wavelengths of interest. Table 1 and Fig. 5 demonstrate that optical crosstalk increases from 
470 nm to 600 nm, because absorption decreases with increasing wavelength, and thus 
photons at longer wavelengths travel further in the tissue. Overall, higher crosstalk is expected 
in normal tissue, since absorption is lower than in malignant tissue. 
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Fig. 5. Contour lines of simulated crosstalk as a function of tissue optical properties for the 
probe geometry reported herein. This crosstalk was simulated using a forward Monte-Carlo 
model for the center pixel of the imaging array illuminated using (a) a reflective tube and (b) an 
absorbing tube. 

The crosstalk in the system was experimentally determined by measuring the diffuse 
reflectance from the two tissue phantoms listed in Table 1 while using beam blocking material 
to block illumination for all PD apertures except the PD pixel under consideration, i.e., only 
one PD aperture was illuminated, and the photocurrent from all PD pixels was measured. 
Beam blocking material (Thorlabs TB4 sliced to 0.6 mm thickness) with a single aperture 



#176730 - $15.00 USD Received 27 Sep 2012; revised 7 Nov 2012; accepted 7 Nov 2012; published 8 Nov 2012 
(C)2012OSA 1 December 2012 /Vol. 3, No. 12 / BIOMEDICAL OPTICS EXPRESS 3219 



allowed illumination of only the PD pixel under consideration. This provided a measure of the 
diffuse reflectance signal at the illuminated PD pixel in the absence of any photons from 
neighboring PD apertures. A second measurement was made by removing the beam block to 
illuminate all PD apertures simultaneously, while keeping the rest of the optical set-up 
unchanged. The two successive measurements were used to calculate the crosstalk using the 
following equation: 



Crosstalk(%) = 



Iph, API — Iph, SPI 
Iph, API 



xlOO, 



where, I p i,_api is the measured PD photocurrent with all photodetector apertures illuminated, 
and I p h,spi is the measured PD photocurrent with a single photodetector pixel illuminated. 

Figure 6(a) and Fig. 6(b) show the worst case scenarios: the center and side pixels with a 
reflective tube (i.e., the tube that results in reflected rays forming large angles with the axis of 
the tube, to emerge from PD apertures) and the absorbing tube for normal, low absorption 
tissue. Although the absorbing tube results in lower optical throughput, with longer ITIA 
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Fig. 6. (a) Crosstalk measured for a center pixel of the imaging array on normal tissue 
mimicking phantom using a reflective and an absorbing tube; (b) Crosstalk measured for a side 
pixel of the imaging array on normal tissue mimicking phantom using a reflective and an 
absorbing tube; (c) Measured and simulated crosstalk for a center pixel of the imaging array on 
a normal tissue mimicking phantom using an absorbing tube; (d) Measured and simulated 
crosstalk for a side pixel of the imaging array on a normal tissue mimicking phantom using an 
absorbing tube. 
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integration time, the SNR and optical crosstalk are superior for the absorbing tube. Figure 6(c) 
and Fig. 6(d) show a comparison between the theoretically predicted crosstalk and 
experimentally observed crosstalk for the absorbing tube (longer integration time) for the 
center (c) and side (d) pixels shown in Fig. 6 (a) and Fig. 6(b), respectively. Two full data sets 
were measured for the center pixel, which varied by 6%-7% in magnitude, but followed the 
same trend (i.e., higher crosstalk with reflective tube, increasing crosstalk at longer 
wavelengths). The data set with best theoretical agreement is shown in Fig. 6(a) and Fig. 6(c). 
The good agreement between the theoretical and measured crosstalk validates the model, and 
offers the potential for this model to be used in future simulations to predict the crosstalk for 
imagers with different pixel configurations. 

5. Comparison to previously developed fiber based DRS imaging systems 

Multipixel fiber-based illumination and collection is a faster method than single probe 
scanning, but poses a challenge in scaling imaging probes for large-area, high-resolution DRS 
imaging. Here, we have reported the first contact-based multipixel DRS imaging system with 
fiberless illumination of tissue, as well as one in which the size of the detector and aperture as 
well as the spacing between adjacent channels can be customized. This system demonstrates a 
nearly 100-fold reduction in size over a previously reported clinically tested bench-top DRS 
imaging system for breast tumor margin assessment, making it amenable for tumor margin 
assessment within the space constraints of an operating room. The resolution of the custom 
imaging array is 4.5 mm compared to 10 mm in the clinically tested multipixel fiber probe 
[12]. In addition, replacement of the imaging spectrograph (-$12,000) and cooled CCD array 
(-$20,000) [13], with eight discrete wavelength filters (-$1440), and a custom imaging array 
(-$1400) with custom read-out electronics (-$1000), respectively, has led to an 8-fold 
reduction in cost compared to the clinically tested bench-top system. Table 2 shows a system 
level comparison of the reported system to the previously reported bench-top system [12]. 



Table 2. System level comparison of the clinically tested bench-top system [12] to the 
system reported herein 
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6. Conclusion 



For the system reported herein, individual system components for light detection (PD array), 
light delivery (free-space light delivery tube), and electronic read-out (8-channel ITIA circuit) 
have been designed, independently optimized, fabricated, and characterized, followed by 
quantification of system SNR and crosstalk using breast tissue mimicking phantoms. The 
system SNR was measured to be 35 dB to 45 dB for X = 470 nm to 600 nm on benign and 
cancerous breast tissue mimicking phantoms using two free-space light delivery strategies; a 
reflective tube that optimized power through PD apertures at the cost of a larger angular 
distribution of rays, and an absorbing tube that results in lower power throughput but also a 
smaller angular distribution of rays through the PD apertures. The angular distribution of rays 
through the PD apertures has implications on the optical crosstalk of the system. Optical 
crosstalk in tissue was measured on breast tissue mimicking phantoms and was compared for 
two different light delivery strategies. Forward Monte-Carlo modeling was conducted to 
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provide a preliminary estimate of system crosstalk, and this model was validated against the 
measured crosstalk. Longer ITIA integration times (100 ms) and use of an absorbing tube for 
light delivery resulted in high SNR (36 dB to 48 dB), while minimizing crosstalk (lower than 
15%) for an imaging resolution of 4.5 mm. 
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